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Using Chandrasekhar's S-type coupling between rotational and intrinsic vortical modes one may 
simply reproduce the HFB dynamical properties of rotating nuclei within Routhian HF calculations 
free of pairing correlations yet constrained on the relevant so-called Kelvin circulation operator. 
From the analogy between magnetic and rotating systems, one derives a model for the quenching of 
pairing correlations with rotation, introducing a critical angular velocity — analogous to the critical 
field in supraconductors — above which pairing vanishes. Taking stock of this usual model, it is 
then shown that the characteristic behavior of the vortical mode angular velocity as a function 
of the global rotation angular velocity can be modelised by a simple two parameter formula, both 
parameters being completely determined from properties of the band-head (zero-spin) HFB solution. 
From calculation in five nuclei, the validity of this modelised Routhian approach is assessed. It is 
clearly shown to be very good in cases where the evolution of rotational properties is only governed 
by the coupling between the global rotation and the pairing-induced intrinsic vortical currents. It 
therefore provides a sound ground base for evaluating the importance of coupling of rotation with 
other modes (shape distortions, quasiparticle degrees of freedom). 



PACS numbers: 21.60.Jz, 21.60.Ev, 21.10.Re 
I. INTRODUCTION 



In a previous paper (hereafter referred as I) we have 
provided some evidence that the dynamical effects of 
pairing correlations at finite spins could be very well 
represented indeed by an intrinsic flow being both non- 
deforming and counter-rotating (with respect to the 
global rotation) pj. This has been borne out from the 
study of rotational bands in three heavy nuclei chosen 
to represent rather different cases (as far as nucleon 
numbers, deformation, spin values or pairing correla- 
tion content are concerned). To achieve this, we have 
performed two sets of microscopic calculations under a 
Routhian-type constraint (namely including in the vari- 
ational quantity for the rotating case a — Q ■ J term with 
usual notation). The first calculations have included 
pairing correlations within the Hartree-Fock-Bogoliubov 
(HFB) formalism. The second ones were of the Hartree- 
Fock type (i.e., without any pairing correlations) with a 
double Routhian-type constraint, as presented in Refs. 
0, 0] and recalled here in Appendix [U] which we will 
refer to below as HF+V calculations . The "measuring 
stick" for the latter was the so-called Kelvin circulation 



operator K x (see, e.g., Ref. 5J) which is well suited for 
the description of collective modes dubbed after Chan- 
drasekhar as S- type ellipsoids . Indeed, upon imposing 
via the Routhian double-constraint that the second type 
of solutions should have the same (Ki) expectation value 
at a given value of the spin I as those corresponding to 
the HFB calculations, it has been demonstrated in I that 
both calculations yielded the same rotational properties. 

As shown also in I, the intensity of such counter- 
rotating currents (measured, e.g., by (Ki)) behaves as 
a function of the global rotation angular velocity in a 
very characteristic way. One may figure out that such a 
reactive mode of the fluid should be proportional to the 
excitation velocity field intensity (i.e., fi). Moreover it 
is also safe to assume that K^i)) should be an increas- 
ing function of the pairing correlations measured by some 
power of —E COII where -E C orr is the (negative) pair corre- 
lation energy. Now, it is well known (see Ref. 0) that 
pairing correlations tendto decrease upon increasing Q. 
When f2 gets larger therefore, one should expect a bal- 
ance between two competing phenomena resulting in a 
maximum of |(lCi)| as a function of f2 between two lim- 
iting cases. One is ft — and the other is obtained for a 
critical value f2 c corresponding to the alleged transition 
between a "superfluid" and a "normal" phase of nuclear 
matter. 
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A word of caution is worth adding at this point. It 



2 



is absolutely not clear that a transition between nor- 
mal and supcrfluid nuclear phases should occur within 
a physically acceptable range of fl values for such finite 
Coulombian systems as atomic nuclei. In this respect it 
is well known that the usual Bogoliubov treatment (or 
its Bogoliubov- Valatin or BCS reduction) tends to over- 
estimate the instability of the superfluid solution in low 
pairing regimes in general, and at high fl values in par- 
ticular. Only approaches which would conserve explicitly 
the particle number, could cast a priori some light on the 
existence of such a transition. Such calculations are cur- 
rently achieved 8] within an extension of the approach 
developed in Ref. 0] to describe correlated states in even- 
even nuclei (there, of course, the time-reversal symmetry 
was not broken). In the present paper however, we will 
overlook this problem and stick to the usual HFB ap- 
proach which entails naturally that this particular part 
of our description in the limit of high f2 regime (i.e., close 
to the critical £l c value) should be considered as rather 
schematic indeed. 
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FIG. 1: Kelvin circulation mean value in ft units for neutrons 
(upper panel) and protons (lower panel) in the three rota- 
tional bands studied in I. The convention in use is the fol- 
lowing: circles for 150 Gd, squares for 192 Hg and triangles for 
254 No. Calculations are performed within the HFB formalism 
for 254 No and within the HFB+LN formalism for 150 Gd and 



The above described competition pattern has been il- 
lustrated in I by curves (see Fig. 1 there) exhibiting (Ki) 
as a function of f2. The main features of these results are 
sketched here in Fig. ^for the isoscalar values of (K\) 
in the three considered rotational bands (yrast superde- 
formed bands of 150 Gd and 192 Hg, and ground-state band 
of 254 No). As seen on this Figure, it is found only in 
the case of the 192 Hg nucleus that the range of VL values 
spanned by our calculations (or more precisely by corre- 
sponding available experimental data) allows to exhibit 
a clean cut extremum in the corresponding curve. For 
the 150 Gd nucleus one merely displays the up-going part 



of the curve (i.e., after the minimum of (K\) and before 
f2 c ) while, quite on the contrary, for the 254 No nucleus, 
one only sees in Fig. ^ the down-going part of the curve 
(i.e., just above £7 = 0) reaching its minimum, only for 
protons. 

In addition, in the 150 Gd case, a very interesting phe- 
nomenon is present and is worth discussing. As presented 
in I, one sees that substantial (K\) values are present for 
both neutrons and protons (namely about +1 h and —2ft. 
respectively) even when pairing correlations have disap- 
peared around KI = 0.6 MeV. This phenomenon, also 
present for the proton Kelvin circulation in 192 Hg (whose 
pairing correlations vanish around KI — 0.4 MeV), is 
clearly a consequence of shell effects as it has been shown 
in I to be also present within purely rotating Routhian 
HF calculations. The behavior of (K\) in No is also 
worth a closer look. Indeed, as the neutron pairing en- 
ergy is almost vanishing (see Fig. 2 in I) even for KI = 0, 
a large negative Kelvin circulation mean value is observed 
in Fig.^for neutrons. This phenomenon also could be at- 
tributed partly to shell effects, but it is probably mostly 
related to a dragging effect of the protons. 

Another way of representing the same trends is to plot, 
as done here in Fig.|3for the three considered nuclei, the 
values of the Lagrange multiplier uj(£1) of the Routhian 
(|C1|) determined in I to obtain the relevant HFB (Ki) 
values as a function of f2. As seen here they exhibit a 
parabolic-like behavior. Being related explicitly to the 
extra-contribution to {K\) due to pairing correlations, it 
is not surprising that such u> (Q) curves are all located in 
our case in the uj < part of the plane. 
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FIG. 2: Kelvin circulation velocity u (in keV) as a function 
of the rotation angular velocity Q in MeV for the three nuclei 
considered in I and calculated within the constrained HF+V 
formalism. The symbol conventions are the same as those in 

Fig.m 

It is the aim of the present paper to understand quan- 
titatively the pattern of these lu(£1) curves. To put it in 
more operational terms we want to find whether we are 
able to deduce from the geometrical and pairing correla- 
tion properties of HFB solutions at zero-spin, the prop- 
erties of rotational states. This of course would be only 
possible in cases where one deals with a "good" rotational 
band behavior, namely wherever the rotation does not 
couple significantly with other degrees of freedom like vi- 
brations (e.g., no rotational stretching or anti-stretching) 
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or single-particle excitations (e.g., no back-bending). 

The paper will be organized as follows. In SectionHTIwc 
will briefly develop a simple magnetic-like model inspired 
from Ref. using the well-known supraconductor criti- 
cal field concept (see, e.g., |l0|) to describe the reduction 
of pairing correlations upon increasing f2. Section lTTTl will 
be devoted to the proposition and discussion of an ansatz 
for the relation between u> and f2. We will present in Sec- 
tion llVI the prescriptions to determine the parameters of 
the above relation from the zero-spin HFB solution of 
a given nucleus. Finally, the validity of the model so 
established will be tested against actual HFB Routhian 
calculations for five nuclei (including the three already 
considered in I) and some conclusions will be drawn in 
Section El 



II. A SIMPLE MAGNETIC-LIKE ANALOG 

In the simple model presented here (which is directly 
inspired from the rotational quenching mechanism of 
pairing correlations proposed by Mottelson and Valatin 
0, the so-called Coriolis anti-pairing effect), we assume 
that the action of a global rotation on a pair of particles 
moving in opposite directions on a given orbital may be 
described roughly in terms of two interacting magnetic 
dipoles (Di, D 2 ) plunged in an external magnetic field 
B. The total energy of such a system is made of a term 
E COII describing the dipole-dipole coupling and a term 
E m taking into account the coupling with the external 
field. However, while the latter is understood as a mag 
netic field-dipole interaction term (see, e.g., Ref. 
Chapter 5) 



-13B-{D 1 +D 2 ) 



(3>0, 



(1) 



the former does not originate, in our case, from a mag- 
netic dipole-dipole interaction (which actually would 
tend to align D\ and D 2 ), but rather from the strong 
interaction favoring quite on the contrary an anti- 
alignment of D\ and Therefore i? cor r (which will 
be later called the correlation energy) is not given by 
the standard magnetic dipole-dipole expression (see, e.g., 
Ref. [Tl|, Chapter 4) but rather, assuming the simplest 
possible form compatible with the physical requirements, 
as: 



lAli^i) 



a > 0. 



(2) 



Assuming B (whose norm will be denoted as B) to lie 
on the 3-axis and Di to belong to the (1,3) plane, one 
gets for the total energy the following expression (where 
di , 4>i are the usual spherical coordinates determining the 
Di direction): 



'tot 



= ad 2 [ sin(#i) sin(02) cos{fa) + cos(0i) cos(02) — l] 
- /?.Bd[ cos(0i) + cos(0 2 )] 

(3) 



where d is the common norm of Di and D 2 . 

The variation of E to t with respect to fa leads to an 
extremum either at sin fa = or else at sin0i and/or 
sin 9 2 = 0. In the latter case, it is easy to check that to 
get an extremal energy, the vanishing of sin 9\ entails the 
vanishing of sin 6 2 and conversely. One is thus left with 
sin 6\ = sin 62 = 0. 

Discarding for a while this solution let us concentrate 
now on the former case. Clearly to get a minimal energy 
fa = 7r is to be preferred to fa = 0. The variation of 
Etot with respect to 0i and 62 leads now to 

/IB sin 0! = f3B sin 2 = ad sin(0i + 2 ) (4) 

which admits two classes of solutions: 



n = t?2 = a , 
h + 6 2 = 7T . 



(5a) 
(5b) 



In the first case, the remaining variational equation 

£Bsin0 = adsin(20) (6) 
has two solutions. The first one is defined by 

(7) 

whose corresponding equilibrium energy is 

I3 2 B 2 



(3B 
2ad 



E u 



-ad- 2 



2a 



(8) 



The second solution is defined by sin0 = 0. Therefore, 
all variational solutions other than the one of eq. Q 
correspond to sin0i = sin 02 = 0, namely : 



VI — u 2 
01 = 02 

h = 0, n 

h = 7T, 



: 
7T 

2 = 7T 
' 2 =0 





= -2(3Bd, 


(9a) 


^tot 


= 2pBd, 


(9b) 


Etot 


= -2ad 2 . 


(9c) 



While searching for the most stable equilibrium solution, 
it is found upon comparing the total energies given in 
eqs. (JSJl and (|9al9cl) that the former is to be retained 
(since in particular ad 2 + (3 2 B 2 / (2a) — 2f3Bd is clearly 
positive) . 

For the solution of eq. (JSJ, as expected, there is an 
upper limit for the norm B c of B given by 



2ad 



(10) 



corresponding to a full alignment (0 = 0). The corre- 
sponding correlation energy E corr at equilibrium assumes 
the following form: 



2ad 2 



(11) 
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One may note, parenthetically, that the critical field 
value is obtained upon equating twice the correlation en- 
ergy without field with the dipole magnetic interaction 
in the aligned case 

-Aad 2 = -f3Bd (12) 

which is very similar in spirit with the way in which one 
determines the critical magnetic field in the classical the- 
ory of supraconductivity (see, e.g., Ref. [lCj)- 

III. MODEL ASSUMPTIONS IN THE CASE OF 
NUCLEAR ROTATIONAL STATES 



For 7^1, one gets two extrema for u> corresponding 
to: 

ni= 7 =%= ; n 2 = n c (16) 

V 1 + 27 

(for 7 = 1 one would have only the variational solution 
corresponding to Oi ) . The correlation energy associated 
with the f2i solution, expressed in units of its value at 
zero angular momentum (f2 = 0) is: 

flcorr(ftl) = 2 7 

E COII (0) l + 2 7 - 1 ' 



As well known, the equations of motion of a particle 
coupled to a magnetic field by its charge, are similar to 
those of a massive particle moving in a rotating (non- 
inertial) frame upon identifying, with a proportionality 
constant, the magnetic field B and the angular velocity 
57. It allows therefore to schematically model the cou- 
pling of a single-particle orbital motion in the rotating 
frame with the global rotation (in the laboratory frame) 
as the magnetic dipole interaction considered in the pre- 
vious section [see eq. |JT]|]. Similarly, the corresponding 
interaction energy between two particles moving in op- 
posite directions on a same orbit could be approximated 
as: 



E r 



X 



^ -1 



(13) 



introducing thus a critical angular velocity J7 C which will 
be defined below. 

The preceding could pertain, a priori, to the descrip- 
tion of the internal and external energies of a single pair. 
It is our contention, however, that such a quantity be- 
haves, up to a multiplicative constant, as the energies of 
all the pairs to be considered in the nucleus. Therefore, 
after suitably redefining the phenomenological parame- 
ter A, one may consider the expression above given in 
eq. (|13l) as the total nuclear pairing correlation energy. 

Now, we are proceeding with the major assumption 
behind our model. Namely, we assume that the vector 
uj associated with the intrinsic vortical motion is pro- 
portional to the global rotation angular velocity f2. The 
corresponding proportionality factor is negative and we 
assume that its absolute value is an increasing function 
of the total correlation energy of eq. (|T3|> . Specifically, 
we propose that the projection u> of the vector uj on the 
vector 17 is given by: 



-a(E COII )m ; a, 7 >0. 



(14) 



Upon inserting in such an ansatz, the correlation energy 
expression of eq. (T3J, one gets for u> 



k > and 17 G [0, Sl c ] . 

(15) 







uj = — fc!7 


M£)1 







IV. DETERMINATION OF THE MODEL 
PARAMETERS 

Let us summarize the modeling which has been per- 
formed so far. With eq. I|15fl . we have a three-parameter 
expression for the intrinsic vorticity parameter uj in terms 
of the global rotation angular velocity 17, namely: 

i) a scale parameter (17 c ) for the abscissa variable 17, 

ii) a scale parameter (fc) for the studied quantity uj, 
and 

iii) a power parameter (7) expressing the nature of the 
pairing correlation dependence of uj. 

The former parameter 17 c may be fixed a priori for each 
considered nucleus as sketched before. In eq. (|12f) . it is 
given upon equating twice the correlation energy without 
field (thus, here, without rotation), with the equivalent of 
the dipole-magnetic coupling energy in the aligned case, 
namely the rotational energy without pairing. We there- 
fore get 



a. 



'4E corr (fi = 0) 



(18) 



In the above Zrig represents the moment of inertia of an 
unpaired nucleus which could be reasonably well repre- 
sented by the rigid-body moment of inertia associated 
with a one-body reduced local density p(r) obtained, e.g., 
in an Hartree-Fock or Hartree-Fock-Bogoliubov calcula- 
tion as: 



m / p(r) (f — (u ■ r)u)' 



d 3 r, 



(19) 



where u is the unit vector 17/17. However, in doing so 
one would introduce an unwanted contribution of shell 
effects and/or pairing correlations to Ziig through the 
density function p(f). Even though the effect of the lat- 
ter is generally rather small we propose here, to remain 
in the spirit of the model (associated with bulk nuclear 
properties), to rather consider the self-consistent semi- 
classical moment of inertia 3etf (within the Extended 
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Thomas- Fermi scheme) of Ref. [IJ using the same effec- 
tive interaction at the relevant deformation. 

In most current Routhian HFB calculations (includ- 
ing ours 

0) 

using the Skyrme interaction for the "nor- 
mal" mean field, there is an inconsistency in the interac- 
tion used in the Hartree-Fock ( "normal" ) part and in the 
part dealing with pairing correlations. As a result, the 
naturally defined correlation energy, namely the energy 
difference between the correlated and the non-correlated 
solutions, is not relevant. Thus we must replace in eq. 
(|f 8fl . the correlation energy E COIY (fl = 0) by an unequiv- 
ocal quantity. The "abnormal" part of the total HFB 
energy — sometimes called (see, e.g., Ref. Q and Ap- 
pendix [SJ the "pair condensation energy", i? con d — could 
be considered in practice in this case, as a more reliable 
index of the amount of pairing correlations. As demon- 
strated in model calculations performed in the three cases 
considered in I according to the approach of Ref. [jj 
which is free from the usual HFB breaking of the par- 
ticle number symmetry, the pair condensation energy is 
found to be roughly equal to twice the pair correlation 
energy (see appendix 1X1 for details). The equation l(T%|) 
yielding the critical value f2 c becomes thus 



l 2E cond (fl = 0) 
3etf 



(20) 



The scale parameter k can be determined in the follow- 
ing way. If we assume that our relation is of any 
global relevance, it should also be valid in the low angu- 
lar velocity regime. Namely, one should have for Q <C tt c 
whatever the choice of the power 7: 



n 



(21) 



In this adiabatic regime, one obtains an expression for 
the total collective kinetic energy in a HF+V approach 
which is quadratic in Q and 10, as with the notation of 
Ref. [3 



E(Q, w) = i (Am 2 + 2Bujfl + Cfl 2 ) 



(22) 



Upon including the above limiting expression for the ratio 
of uj and f2 one gets the expectation value of the total 
angular momentum as 



(Ji) = 



dE(n,t 



uj— — kQ 



= (C-kB)Q, 



(23) 



which in turn yields the value of the dynamical moment 
of inertia within our model: 



(2) = d(JQ 

Omod dn 



= C-kB 



(24) 



as shown in Appendix IH1 

Now, we define the HFB Routhian dynamical moment 
of inertia 3hfb (in the low angular velocity regime) and 
require it has the same value as Z 



Jhfb 



(Ji) 
(I 



(2) 
mod ' 

? (2) 
J mod ' 



(25) 



One is left therefore, now, with the problem of deter- 
mining the moments A, B and C . A natural way to do 
so is of course to calculate the total energy surface as a 
function of Q and u> in the vicinity of zero for both param- 
eters through doubly constrained HF calculations (called 
in I Routhian HF + V calculations) and then perform 
a quadratic fit. This has been done indeed for the five 
considered cases (see Appendix iBl for the results). Even 
though this does not amount to performing very much 
time-consuming calculations, we have developed another 
approach which is much simpler and seems to be quite 
sufficient for the model calculations tested here. 

One may evaluate from the single-particle eigenstates 
of the corresponding Hartree-Fock solution, the mass pa- 
rameters A, B and C through the Inglis cranking formula 
|l6l Il7j ] where the cranking operators are (see Ref. |l5j) 
the first components of the total angular momentum J\ 
and of the Kelvin circulation operator K\. With usual 
notation one has thus 



C = 2£ 



\(p\Ji\h)\ 

(p — e/i 



B = 2J2 

p . h 



(p — e/i 



(26) 



A = 2J2 



\<p\Ki\h)\< 
£ P — e/i 



/.) . h 

However it is well known (see e.g. Ref. 0|) that 
one should not compare dynamical moments of inertia 
obtained in a Routhian approach with comparable mo- 
ments in an Inglis cranking approach. Indeed the latter 
lack the so-called Thouless-Valatin \lQ terms which are 
merely coming from the time-odd density response to the 
self-consistent time-odd Hartree-Fock mean field. For in- 
stance, it has been estimated in Ref. (2j| for the usual 
cranking moment C and in the presence of pairing cor- 
relations that one gets with the Inglis cranking formula 
an under- estimation of the Routhian value of about |. 

In view of the similarity of the two operators J\ and K\ 
it is reasonable to further assume that the same multi- 
plying factor should be applied also for the mass param- 
eters A and B. These approximations and the value of 
the Thouless-Valatin correction parameter r\ [such that 
the actual values should be enhanced from their Inglis 
cranking crude estimates by a factor (1 + 77)] have been 
assessed from doubly constrained HF calculations in the 
five considered cases in Appendix [Bl As we will see in 
the next section, two different Skyrme effective interac- 
tions have been used in our calculations (SIII [2l| and 
SkM* 22]). The corresponding values for 77 are 0.2 and 
0.1 respectively. 

The difference between these values as well as with 
the value proposed in Ref. [2(j may have two origins. 
The first is that the interactions used in the various self- 
consistent calculations are different. As shown in |23| . 
the Thouless-Valatin corrective terms for Skyrme inter- 
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actions increase when the difference between the real nu- 
cleon mass (to) and its effective value in nuclear matter 
(to*) increases. If this evidence for Skyrme interactions 
is to be extended to the Gogny force in use in [2(j , it ap- 
pears that at least a part of this difference may be due to 
the fact that the ratio to*/to equals 0.70 for the consid- 
ered Gogny force, 0.75 for the SIII force, and 0.79 for the 
SkM* force. The second is that in Ref. [2fJ one has per- 
formed Routhian calculation with pairing correlations. 
As a result one gets the following value for k: 

B 

The last parameter whose value has to be pinpointed 
is the power 7 of the correlation energy dependence of 
u. By definition if 7 = 1 the parameter to depends lin- 
early on the pairing energy (correlation or condensation 
energies) whereas if 7 = 0.5 it would depend linearly on 
the average pairing gap. We have found no binding ar- 
guments to decide a priori on this value. Besides, as seen 
from equation (|17fl in the 7 = 0.5-1 range the location 
of the maximum value of \lu\ is not violently contingent 
upon a specific figure for 7. We have therefore decided 
not to use the freedom of fitting this parameter and there- 
fore we have somewhat arbitrarily decided to take 7 = 1. 

V. RESULTS AND CONCLUSIONS 

We have applied the previously described protocol to 
define the model parameters fl c and k [see eq. <|15[) 
with 7=1] for rotational bands in five nuclei, includ- 
ing the three cases studied in I and adding the ground- 
state bands of 154 Sm and 178 Hf. The Kelvin circulation 
mean values for protons and neutrons calculated in the 
Routhian HFB approach are plotted on Fig. [3] as a func- 
tion of for these two nuclei. The different trends of 
these curves for 154 Sm as opposed to 178 Hf do not reflect 
the behavior (illustrated on Fig. 0} of the condensation 
energies, but some specific dynamical properties which 
will be discussed bellow. 

It is generally considered (see, e.g., 0) that the SIII 
|2lj interaction provides extremely good spectroscopic 
properties for normally deformed nuclei, while the SkM* 
[22| interaction, due to its better surface tension, is well 
suited to the description of superdeformed or heavy (i.e., 
with a fissility close to 1) nuclei. Therefore we have used 
the former for the calculations of the ground-state bands 
of 154 Sm and 178 Hf, and the latter for the three others 
nuclei already considered in I. 

First, we have obtained values of the critical angular 
velocities f2 c according to eq. 120fl . 

For the five nuclei, HFB calculations at zero-spin have 
been performed to yield the corresponding pair conden- 
sation energies -E con d(^ = 0). They have already been 
shown on Fig. 2 in I for 192 Hg and 254 No. In the 
150 Gd case, we have "artificially" constructed a zero- 
spin solution consistent with the deformation of the su- 
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FIG. 3: Same as Fig. [T] for 154 Sm (diamonds) and 178 Hf 
(stars). 
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FIG. 4: Condensation energies (in MeV) as functions of the 
angular velocity Q (in MeV) for 154 Sm (upper panel) and 
178 Hf (lower panel). Proton (neutron resp.) condensation 
energies are represented as a dashed line (dotted line resp.) 
and the total condensation energy is represented as a full line. 
Calculations are performed within the HFB formalism. 

perdeformed band under consideration by constraining 
the HFB solution to have the value of the axial mass 
quadrupole moment of the SD rotational band states (of 
course, without constraint one would have obtained the 
normally deformed 150 Gd equilibrium solution). The re- 
tained values of -E CO nd(^ = 0) are listed in table U 

The values of the semiclassical moment of inertia to 
be used in eq. l{2T)ll (see above the discussion in Section 
IIV|) are also listed in table |U together with the resulting 
il c values. It is encouraging to note that these values are 
actually rather close to the expected values of f2 which 
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FIG. 5: Results of our calculations for 192 Hg (upper panels) and 254 No (lower panels). From left to right we show (as functions 
of the rotation angular velocity Q in MeV) for each nuclei the kinetic moment of inertia and the dynamical moment of inertia 
(in h 2 MeV" 1 units) as well as the Kelvin circulation velocity u> (in keV). The results obtained within various formalisms 
are represented as follows: dotted line for the (pure) HF formalism, full line (dashed line resp.) for the HFB (HF+V resp.) 
formalisms and full line with opened circles for our model. Experimental data for the moments of inertia are represented as 
filled circles. The Lipkin-Nogami correction has been applied for the Hg isotope. 



TABLE I: Input parameter E^ond (in MeV), 3etf and 3hfb 
(in h 2 MeV -1 units) for our model (see text). The output 
parameter Ml c (in MeV) and k are also given. 





•^cond 


3etf 


3hfb 


nn c 


k 


lbu Gd 


10-47 


88-47 


69-74 


0-4864 


0.3443 


154 Sm 


10-03 


77-85 


32-23 


0-5077 


0.6629 


178 ffi 


16-14 


97-77 


26-14 


0-5745 


0.6285 


192 Hg 


6-99 


120-08 


85-94 


0-3411 


0.3159 


254 No 


9-45 


164-53 


68-45 


0-3389 


0.3196 



would correspond to a vanishing of the isoscalar pair cor- 
relation energy as exhibited (see Fig. 2 of I and Fig. |3| 
here) . 

The values of A, B, C and rj have been discussed and 
given in Appendix [5] The values of Zhf b have been 
calculated in the close vicinity of = and may be de- 
duced from the plots of Z^ 2 ' of Fig. 1 in I for 192 Hg and 
254 No. Again in the case of 150 Gd this value is yielded 
by a Routhian HFB calculation with an appropriate con- 
straint on the axial mass quadrupole moment (as above 
discussed). The values of these moments are also listed 
in table [I] together with the resulting values of the k pa- 
rameter. It is worth noting in Table [I] that the k values 
corresponding to a given effective force are not very much 
varying from one nucleus to another. However they seem 
to be strongly dependent of the considered force (about 
0.3 for SkM* and 0.6 for SIII). We do not yet understand 
this feature. 

As a result, the curves uj(SY) yielded by the above de- 



fined parameters are plotted on Figs. EH3 m compari- 
son with the self-consistently calculated ones (HF+V). 
The global agreement is rather satisfactory in a qualita- 
tive fashion for 150 Gd, 192 Hg and 254 No. This is not so 
for 154 Sm and 178 Hf where a sudden raise of u> occurs. 
In these three figures, experimental data are taken from 
Refs. for 154 Sm and 178 Hf, [H for 150 Gd, for 
192 Hg and HHH for 254 No. 

The final step of the assessment of our model consists of 
course in performing Routhian HF+V calculations with 
the model values of u>(f2) and compare their results with 
those of corresponding Routhian HFB calculations. They 
should also be compared with the results of Routhian 
HF+V calculations where the constrained (Ki) expec- 
tation values are those obtained in Routhian HFB cal- 
culations. However we knew already from I (see Fig. 1 
therein) that both are very close indeed. 

The results of the comparison fall into three categories. 
For the 192 Hg and 254 No nuclei (Fig. |3J), the agree- 
ment between the HFB (and thus HF+V) and the model 
results are very good indeed. One exception however 
should be noted, it concerns the high spin part of the su- 
perdeformed band of 192 Hg where pairing correlations are 
strongly damped or disappearing in the HFB approach 
and therefore the simple ansatz here is obviously lacking 
any sound basis. 

The 154 Sm and 178 Hf cases (Fig. |SJ| are interesting in 
that HFB calculations exhibit a up-bending pattern for 
both an d ■ The fact that it corresponds very 
well to the experimental situation in the latter and is 
completely absent in the data for the former is irrelevant 
for our discussion here. It is remarkable that our model 
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FIG. 6: Same as Fig. EJfor 154 Sm (upper panels) and 178 Hf (lower panels). 




0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0., 

HQ. (MeV) HQ (MeV) HQ. (MeV) 



FIG. 7: Same as Fig. E]for 150 Gd. The Lipkin-Nogami correction has been applied. 



just ignores these bendings. Therefore we must conclude 
that these patterns are not to be attributed to a pairing- 
rotation coupling but to something else which is (as well 
known) the intrusion of quasiparticle degree of freedom 
in the rotational dynamics. 

As for the last nucleus of our study ( 150 Gd), Fig. 0dis- 
plays clearly in its 3^ and 3^ parts that obviously the 
quasiparticle degrees of freedom are dominating there in 
explaining the rotational behavior within this superde- 
formed band. Nevertheless, it is to be noted that our 
model works fairly well in the low angular velocity regime 
and only fails above Ml ~ 0.4 MeV after the occurrence 
of a change in the ordering of quasiparticle states. 

Let us summarize what has been learned from such a 
comparison. Actually, one should conclude at two differ- 
ent levels. 

First of all, it appears in some of the studied cases 
that one could skip Routhian HFB calculations at finite 
values of the angular momentum and replace them by 
more handy constrained Routhian HF calculations. For 
that one needs only to perform a HFB calculation in the 
close vicinity of CI = 0. This result is of some practical 
value in that it allows to predict quickly what could be 



the trend of the moments of inertia with respect to the 
angular velocity. 

However, the limitation of such an approach is clear. 
It is only valid in cases where the rotational (global) 
collective mode does not couple either with deformation 
modes (no rotational stretching, anti-stretching or triax- 
ial modes, etc.) nor with single-particle degrees of free- 
dom (no back-bending for instance). It is only able to 
give an account which turns out to be rather good quan- 
titavely of the coupling of the pairing degrees of freedom 
with the rotational mode. 

Taking at face value the rough assumption made in 
previous papers Q, Q that the collective effects of these 
correlations could be mocked up by Chandrasekhar S- 
type ellipsoid velocity fields together with the Mottelson- 
Valatin picture for the rotational quenching of pair- 
ing correlations, our model provides a good reproduction 
of fully self-consistent HFB results. It offers thus rea- 
sonable grounds to assume that both conjectures rather 
accurately describe the microscopic coupling at work in 
rigidly deformed purely collective rotational bands. 
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APPENDIX A: COMPARISON BETWEEN 
CORRELATION AND CONDENSATION 
ENERGIES 

Within the HFB formalism, using different interactions 
in the particle-hole and particle-particle channel, one 
cannot obtain a relevant definition of the pair correlation 
energy which should be the energy difference between the 
correlated (HFB) and uncorrelated (HF) states. On the 
contrary, the pair condensation energy which is defined 
as the "abnormal" part of the HFB energy proportional 
to 7>(kA), where k is the "abnormal" density matrix 
and A is the pairing tensor, is free of any ambiguity. 




j.j - 



O.J 1 1.5 2 

h-intensity 

FIG. 8: Ratio of the condensation energy with respect to the 
correlation energy as a function of the pairing force intensity 
(relative to a somewhat arbitrary nominal value) for the three 
nuclear states considered in I. The symbol conventions are the 
same as those in Fig. Q 



ocal in the case considered bellow where the vacuum is 
not self-consistently defined from the correlated p matrix 
(i.e., the non self-consistent case of Ref. @). 

We have performed HTDA calculations for the three 
nuclear states considered in I, using a <5-type pairing in- 
teraction with intensities varying from 0.2 to 2 times the 
nominal intensity given in Ref. [3(j ■ The ratio of the pair 
condensation energy with respect to the pair correlation 
energy obtained within this formalism plotted in Fig. |S| 
is seen to change rather slightly around 2 as the pairing 
interaction varies. 



APPENDIX B: MASS PARAMETERS WITHIN 
THE HF+V FORMALISM 

For the three considered nuclear states we have com- 
puted the collective kinetic energy surface on a 4x4 mesh 
for angular velocities fl and u> lower in absolute value 
than 30 keV/i -1 . This energy surface has then been fit- 
ted by the quadratic expression of eq. j22jl. The mass 
parameters A, B and C so obtained do take into ac- 
count the Thouless-Valatin self-consistent contributions 
and are listed in table |Hj The Inglis approximation of 
these mass parameters as obtained by eqs. l|2f)|) are also 
given. 



TABLE II: Mass parameters of the collective energy in HF+V 
calculations for the five considered nuclear states in ft 2 MeV -1 
units. 







Inglis 




Thouless-Valatin 




C 


B 


A 


C 


B 


A 


150 Gd 


88-52 


74-79 


86-49 


96-92 


81-40 


94-16 


154 Sm 


67-26 


62-57 


63-47 


84-92 


78-81 


79-14 


178 Hf 


48-17 


43-37 


44-24 


57-40 


51-38 


51-88 


192 Hg 


106-49 


92-26 


101-46 


115-88 


99-27 


108-65 


254 No 


86-70 


78-51 


81-48 


96-95 


87-50 


89-98 



An alternate formalism to describe pairing correlations 
in nuclei, dubbed as "Higher Tamm-DancofT Approx- 
imation" (HTDA), free of particle number symmetry- 
breaking, has been recently presented in Ref. |!j ■ Within 
this formalism, the correlated state is obtained — through 
the diagonalisation of a residual interaction in a many- 
body basis — as a weighted sum of wavefunctions con- 
structed as particle-hole excitations of a Slater determi- 
nant vacuum. There, it is possible to define consistently 
a pair correlation energy as well as a pair condensation 
energy. The former is the difference between the HTDA 
energy and the HF energy of the vacuum Slater determi- 
nant, while the latter is the difference between the HTDA 
energy and the HF energy functional associated to the 
correlated one-body reduced density. Actually, this defi- 
nition of the pair condensation energy is totally unequiv- 



TABLE III: Enhancement of the Thouless-Valatin mass pa- 
rameters relative to the Inglis ones (in percent) within the 



HF+V formalism. 

C B A 

150 Gd 9.5 % 8.8 % 8.9 % 

154 Sm 26.3 % 26.0 % 24.7 % 

178 Hf 19.2 % 18.5 % 17.3 % 

192 Hg 8.8 % 7.6 % 7.1 % 

254 No 11.8 % 11.5 % 10.4 % 



As expected, the Thouless-Valatin values are higher 
than those obtained with the Inglis approximation. The 
enhancement (due to the taking into account of all time- 
odd self-consistent contributions) is shown in table II I II 
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In view of the values exhibited in this table, we can as- 
sume that for the three considered nuclear states, the 
Thouless-Valatin correction roughly enhance the Inglis 
mass parameters by 10 % for the SKM* interaction and 
20 % for SIII. 



APPENDIX C: MOMENT OF INERTIA WITHIN 
DOUBLY-CONSTRAINED ROUTHIAN 
CALCULATIONS 



Within a purely rotational formalism, only the first term 
of the right hand side would be present and would be 
identified with the dynamical moment of inertia. The 
second term which arises here is therefore directly related 
to the added intrinsic vortical mode and one is led to the 
conclusion that the dynamical moment of inertia in our 
model should be taken as 



The doubly-constrained Routhian in use in HF+V 
calculations within the present model writes (see Refs. 
Hill): 



R= H - QJi - u(n)K 1 



(CI) 



The stationary condition S(R) — 0, at a fixed value of f2, 
entails 



l d(g) 



d (Ji) oj(0) d (Kx) 



n dn 



(C2) 



5(2) 



d(Ji) 
dfl 



(C3) 



This is actually what had been assumed in I where it was 
substantiated by the fact that requiring identical values 
of (Jx) and {Ki) within both HFB and HF+V formalisms 
yielded the same angular velocity f2 in the two calcula- 
tions and thus the same moments of inertia. 



[1] H. Laftchiev, D. Samsoen, P. Quentin, and I. N. 

Mikhailov, Phys. Rev. C 67, 014307 (2003). 
[2] I. N. Mikhailov, P. Quentin, and D. Samsoen, Nucl. Phys. 

A627, 259 (1997). 
[3] D. Samsoen, P. Quentin, and J. Bartel, Nucl. Phys. 

A652, 34 (1999). 
[4] D. Samsoen, P. Quentin, and I. N. Mikhailov, Phys. Rev. 

C 60, 014301 (1999). 
[5] G. Rosensteel, Phys. Rev. C 46, 1818 (1992). 
[6] S. Chandrasekhar, Ellipsoidal Figures of Equilibrium 

(Dover, New York, 1987). 
[7] B. R. Mottelson and J. G. Valatin, Phys. Rev. Lett. 5, 

511 (1960). 

[8] H. Lafchiev, J. Libert, and P. Quentin, in preparation. 
[9] N. Pillet, P. Quentin, and J. Libert, Nucl. Phys. A697, 
141 (2002). 

[10] C. Kittel, Quantum Theory of Solids (Wiley, New York, 
1963). 

[11] J. D. Jackson, Classical Electrodynamics (Wiley, New 
York, 1962). 

[12] K. Bencheikh, P. Quentin, and J. Bartel, Nucl. Phys. 

A571, 518 (1994). 
[13] H. Laftchiev, D. Samsoen, P. Quentin, and J. Piperova, 

Eur. Phys. J A 12, 155 (2001). 
[14] H. Flocard, P. Quentin, A. K. Kerman, and D. Vautherin, 

Nucl. Phys. A203, 433 (1973). 



[15] I. N. Mikhailov and P. Quentin, Phys. Rev. Lett. 74, 

3336 (1995). 
[16] D. R. Inglis, Phys. Rev. 96, 1059 (1954). 
[17] D. R. Inglis, Phys. Rev 103, 1786 (1956). 
[18] M. J. Giannoni and P. Quentin, Phys. Rev. C 21, 2060 

(1980). 

[19] D. J. Thouless and J. G. Valatin, Nucl. Phys. 31, 211 
(1962). 

[20] J. Libert, M. Girod, and J.-P. Delaroche, Phys. Rev. C 

60, 054301 (1999). 
[21] M. Beiner, H. Flocard, N. Van Giai, and P. Quentin, 

Nucl. Phys. A238, 29 (1975). 
[22] J. Bartel, P. Quentin, M. Brack, C. Guet, and H.-B. 

Hakansson, Nucl. Phys. A386, 79 (1982). 
[23] M. J. Giannoni and P. Quentin, Phys. Rev. C 21, 2076 

(1980). 

[24] R. B. Firestone et al., Table of isotopes, eighth edition 
(J. Wiley and sons, New York, 1996), and Refs. quotted 
therein. 

[25] P. Fallon et al., Phys. Lett. 257B, 269 (1991). 
[26] B. Gall et al., Z. Phys. A 347, 223 (1994). 
[27] P. Reiter et al., Phys. Rev. Lett. 82, 509 (1999). 
[28] M. Leino et al., Eur. Phys. J. A 6, 63 (1999). 
[29] P. Reiter et al., Phys. Rev. Lett. 84, 3542 (2000). 
[30] T. L. Ha, P. Quentin, and D. D. Strottman, in prepara- 
tion. 



